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Nodal noncentrosymmetric superconductors have topologically nontrivial properties manifested by protected 
zero-energy surface states. Specifically, it was recently found that zero-energy surface flat bands of topological 
origin appear at their surface. We show that the presence of certain inversion-type lattice symmetries can give 
rise to additional topological features of the gap nodes, resulting in surface states forming one-dimensional 
arcs connecting the projections of two nodal rings. In addition, we demonstrate that Majorana surface states 
can appear at time-reversal-invariant momenta of the surface Brillouin zone, even when the system is not fully 
gapped in the bulk. Within a continuum theory we derive the topological invariants that protect these different 
types of zero-energy surface states. We independently derive general conditions for the existence of zero-energy 
surface bound states using the complementary quasiclassical scattering theory, explicitly taking into account 
the effects of spin-orbit splitting of the bands. We compute surface bound-state spectra for various crystal 
point-group symmetries and orbital-angular-momentum pairing states. Finally, we examine the signatures of 
the arc surface states and of the zero-energy surface flat bands in tunneling-conductance spectra and dicuss how 
topological phase transitions in noncentrosymmetric superconductors could be observed in experiments. 
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I. INTRODUCTION 

Systems with strong spin-orbit coupling (SOC) have re- 
cently attracted a great deal of attention. Prime examples 
are topological insulators, where strong spin-orbit interactions 
give rise to a non-trivial band topology, leading to topolog- 
ically protected zero-energy surface states.GB Another class 
of compounds for which SOC plays an important role are 
superconductors without inversion symmetry.'*^ In these re- 
markable materials, Rashba-type antisymmetric spin-orbit in- 
teractions lift the spin degeneracy of the electronic bands and 
generate complex spin textures in the electron Bloch func- 
tions. In the superconducting state, the antisymmetric SOC 
gives rise to the admixture of even-parity spin-singlet and 
odd-parity spin-triplet pairing components and, importantly, 
allows a non-trivial topology of the Bogoliubov-quasiparticle 
wavefunctions.^ ^ Akin to topological insulators, this non- 
trivial wavefunction topology results in various types of pro- 
tected zero-energy states at the edge or surface of noncen- 
trosymmetric superconductors (NCSs).^ For instance, a 
fully gapped NCS with nontrivial topology supports linearly 
dispersing helical Majorana modes at its boundary.**^ In 
three-dimensional systems, the stability of these Majorana 
surface states is protected by an integer (Z) topological invari- 
ant, i.e., the three-dimensional winding number,'' whereas 
in two-dimensional systems a binary (S2) to pologic al number 
guarantees the robustness of the edge modes.'^^^EBES 

Remarkably, topologically protected zero-energy boundary 
modes also occur in NCSs with line nodes.l^tSl in particu- 
lar, it has recently been shown that dispersionless zero-energy 
states (i.e., flat bands) of topological origin generically ap- 
pear at the surface of three-dimensional nodal NCSs.l2E2lThese 
zero-energy flat bands are confined to regions of the two- 
dimensional surface Brillouin zone (BZ) that are bounded 
by the projections of the nodal lines of the bulk gap.^"^ The 
topological protection of these dispersionless boundary states 



is linked to the topological characteristics of the nodal gap 
structure via a bulk-boundary correspondence. In fact, the 
stability of both the dispersionless zero-energy surface states 
and the line nodes of the bulk gap is ensured by the con- 
servation of the same integer topological invariant, namely 
the one-dimensional winding number Apart from these two- 
dimensional surface flat bands, certain NCSs also support 
zero-energy boundary states that form one-dimensional open 
arcs i n the su rface BZ, connecting the projection of two nodal 
rings. "-^1221231111 Moreover, it has recently been reported that 
Majorana surface states ca n oc cur at time-reversal-invariant 
momenta of the surface BZ,'^G3 even if the superconductor is 
not fully gapped in the bulk. 

Topological surface states are generic features of NCSs 
whose spin-triplet pairing component is at least as strong 
as the spin-singlet one. It is therefore quite reasonable 
to expect that these zero-energy boundary states occur in 
the superconducting state of the heavy fermion compounds 
CePt3Si,^CeRhSi3,EIIand CeIrSi3,'^as weU as in YsCg,^! 
Li2Pd^Pt3_^B,l^3® and M03 AlaCPH All of these noncen- 
trosymmetric materials show strong spin-orbit interactions, 
with the spin-orbit band splitting far exceeding the supercon- 
ducting energy scale.l^lFor the study of surface phenomena in 
NCSs it is thus important to explicitly account for the strong 
SOCl'^'^^'^^EEl a subtlety that has been overlooked in many 
previous works .G^EHSlli 

The main aim of this paper is to provide a comprehensive 
classification of the topological features of three-dimensional 
nodal NCSs. In particular, we demonstrate that the presence 
of certain inversion-type symmetries gives rise to additional 
topological features of the gap line nodes. These features are 
characterized by a two-dimensional topological invariant 
and manifest themselves on the surface as one-dimensional 
arcs of zero-energy states, terminating at the projection of the 
nodal lines onto the surface BZ. We also examine the appear- 
ance of Majorana surface states at time-reversal-invariant mo- 
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menta of the surface BZ and show that the topological prop- 
erties of these linearly dispersing modes are described by a 
one-dimensional invariant. Using a continuum model of 
NCSs, we derive expressions for the invariants and the 
winding number that protect the Majorana modes, the arc sur- 
face states, and the surface flat bands, respectively. 

We illustrate these topological features by investigating the 
surface-bound-state spectrum and the tunneling conductance 
of NCSs using quasiclassical scattering theory. Allowing for 
non-negligible SOC, we consider three different experimen- 
tally relevant crystal point-group symmetries and explore the 
effects of pairing with higher orbital angular momenta [e.g., 
(d^2_y2 + /)-wave pairing]. We find that higher-orbital-an- 
gular-momentum pairing leads to additional topologically sta- 
ble line nodes in the bulk gap. Correspondingly, there appear 
extra surface flat bands associated with these additional nodal 
lines. We show that the surface flat bands and the arc states 
leave unique signatures in the tunneling-conductance spec- 
tra. Finally, we investigate topological phase transitions in 
NCSs, i.e., zero-temperature quantum phase transitions where 
the momentum-space topology of the quasiparticle spectrum 
changes abruptly as a function of the singlet-to-triplet ratio in 
the pairing amplitude. We argue that anomalies in the density 
of states and the tunneling conductance provide experimental 
fingerprints of these zero- temperature phase transitions. 

The remainder of this paper is organized as follows. Sec- 
tion |ll] discusses the model Hamiltonian and its symmetries. 
In Sec. Ill we derive the topological invariants characterizing 
topological properties of both the nodal lines and the surface 
states and give a detailed discussion of the topological crite- 
ria for the existence of zero-energy surface states. By use of 



quasiclassical scattering theory, we derive in Sec. IV general 
conditions for the occurrence of surface bound states in terms 
of sign changes of the superconducting gap functions across 
the Fermi surface. We show that these conditions are in per- 



fect agreement with the topological criteria given in Sec. Ill 
In addition, we present in Sec . |IV | surface-bound state spectra 
for three different crystal point-group symmetries and for var- 



ious surface orientations. In Sec. V A we compute the tunnel- 
ing conductance between a normal metal and an NCS, identify 
the signatures of the zero-energy surface flat bands and the arc 
surface states in the tunneling spectra, and discuss topological 
phase transitions. Our conclusions and outlook are given in 
Sec.NB 



II. MODEL HAMILTONIAN AND SYMMETRIES 

We consider a three-dimensional single-band BCS su- 
perconductor with noncentrosymmetric crystal structure and 
Rashba-type SOC. On a phenomenological level, such a 
superconductor is described by the Bogoliubov-de Gennes 
Hamiltonian H = ^J2k *k'^(k)*k' with 



_ / h{k) A(k) 



(la) 



and *k = (Ckt'Ck;,cLk^,c^-ki)'^' where 4^ (c^^) denotes 
the electron creation (annihilation) operator with momentum 



k and spin a. The normal-state dispersion of the electrons in 
the spin basis is given by 



/i(k) = ekO-Q + gk • cr, 



(lb) 



where £k = h^i^^/{2m) — /i, gk denotes the SOC potential, 
(7 ~ {ax, (Jy, Cz)^ are the three Pauli matrices, and cto stands 
for the 2x2 unit matrix. In the so-called helicity basis the 
normal-state Hamiltonian ( [Tb] l takes diagonal form, /i(k) = 
diag(Ck 7 Ck where = Sk ± |gk| is the dispersion of the 
positive-helicity and negative-helicity bands, respectively. 

Due to the lack of inversion symmetry, the superconduct- 
ing gap generally contains an admixture of even-parity spin- 
singlet and odd-parity spin-triplet pairing states. 



A(k) = (V'kcro + dk • cr) {i(7y) , 



(Ic) 



where ip]^ and dk represent the spin-singlet and spin-triplet 
components, respectively. It is well known that in the absence 
of interband pairing, the superconducting transition tempera- 
ture is maximized when the spin-triplet pairing vector dk is 
aligned with the polarization vector gk of the SOC*^ Hence, 
we parametrize the singlet and triplet components of the su- 
perconducting gap function as 



i^k - A,/(k) = 



Ao/(k), 



r + l 

dk = At/(k) lk = ^Ao/(k) Ik, 

r + l 



(2a) 



(2b) 



where Ik = gk/A, with A the SOC strength. Here, r = 
As/ At denotes the ratio between the singlet and triplet pair- 
ing components. The pairing amplitudes A^ and A( are as- 
sumed to be positive and constant, i.e., r > 0. We have 
included in Eq. (|2]i a structure factor /(k), which allows us 
to investigate the effects of higher-orbital-angular-momentum 
pairing. In the following, we consider the three cases 



/(k) 



for (s + p)-wave, 

for {dx2^y2 + /)-wave, (3) 

for (dxy + p)-wave, 



where kp is the Fermi vector in the absence of spin-orbit in- 
teractions. With Eq. it follows that the gaps on the two 
helicity bands are given by A^ = /(k)Ao(?' ± |lk|)/(r + 1). 

The specific form of the SOC vector gk (and hence of Ik 
and dk) is constrained by time-reversal symmetry and the 
crystallographic point-group symmetries of the superconduc- 
tor Time-reversal symmetry requires gk to be real and an 
odd function of k. An element R of the crystallographic point 
group Q of the NCS can be represented as either a proper (for 
det R = +1) or an improper (for det i? = — 1) rotation. Thus, 
the Bogoliubov-de Gennes Hamiltonian ([T]) transforms under 
an operation i? of ^ as 



Uj,n{R-'-k)ul 



-Hik) 



(4) 



with R = det{R)R, [/^ = diag(uj^, u*^), and the spinor 
representation of R, i.e., = exp [—i{d/2)h ■ cr]. Here, fi 
denotes the unit vector along the rotation axis of R and 6 is 
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the angle of rotation. It follows from Eq. ([T]i that the lattice 
symmetries Q impose the constraint 



gk = det{R)Rgji-ii,. 



(5) 



on gk. To determine the form of gk, we employ a small- 
momentum expansion. Using Eq. (|5]l we find that for the 
tetragonal point group Q = (relevant for CePtaSi, 

CeRhSis, and CelrSia) the lowest-order symmetry-allowed 
term is^l 



gk = A(fcyX - k^y) 



(6a) 



For the cubic point group Q = O (represented by 
Li2Pdj;Pt3_j;B and M03AI2C) the vector gk takes the form 



gk = A 



fc,(l+g2 [kl + kl])x + ky{l + g2 [kl + kl])y 
k,{l+g2[kl + kl])z\ (6b) 



with the second-order SOC g2- Finally, for the tetrahedral 
point group Q = Td (relevant for Y2C3) the small-momentum 
expansion of gk yields 



gk 



A 



k,{kl-kl)^+ky{kl~kl)y + kM-kl)i 



(6c) 



Before deriving the relevant topological invariants for 
Hamiltonian ([TJ, we first discuss in some detail the dis- 
crete symmetries responsible for the nontrivial topological 
characteristics of 'H(k). According to the classification of 
Refs. T5I, I^T, and"??, 'H(k) belongs to symmetry class Dili 
since it satisfies two independent antiunitary discrete sym- 
metries: Particle-hole symmetry (PHS) C = ICUc, with 
Uc — (7i CTo,^^ and time-reversal symmetry (TRS) T = 
ICUt, with Ut = ctq (8) icr2, where K, denotes the complex 
conjugation operator TRS acts on the Bogoliubov-de Gennes 
Hamiltonian T-L{k) as 



UTn^{-'k)ul = +n{k), 



and PHS as 



UcH'{~k)U^ 



c 



(7) 



(8) 



Combining TRS and PHS, one obtains a third discrete sym- 
metry, the so-called chiral symmetry, which acts as 



uln{k)Us = -n{k) 



(9) 



with the unitary matrix Us — iUtUc — —(^i ® ^2- As 
we will explain in Sec. Ill it is the chiral symmetry (|9| that 
leads to the protection of the zero-energy surface flat bands. 
Since 'H(k) anticommutes with the unitary matrix Us, it can 
be brought into block off-diagonal form. This is achieved by a 
unitary transformation that diagonalizes Us, e.g., WUsW'^ = 
diag(CTo, —<^o), with 



W = 



1 



-i(72 



The transformed Hamiltonian reads 



-0-2 

-zero 



D, 

dI 



(10a) 



(10b) 



where the off-diagonal block is given by 

-Dk = (BkO-Q + v4k Ik ■ cr) {-1(72) 



(10c) 



with the short-hand notation Ak = A + iAt/(k) and Bk = 
Ek + *As/(k). We will see in the next section that the topo- 
logical invariants characterizing the topology properties of 
'H(k) can be conveniently defined in terms of the off-diagonal 
block Dk, Eq. ( [TOj l, or its flat-band version. It is found that 
a crystallographic point-group operation R acts on Dk as 
upjDji-i]^u^p, ~ Z?k, whereas TRS implies Dk = ^^-k- 

Depending on the point group Q and the specific form of 
gk, the Hamiltonian Hik) may in addition to Eq. ^ also sat- 
isfy another "time-reversal"-like symmetry that acts only on a 
two-dimensional plane within the three-dimensional BZ, i.e. 
a symmetry that acts on Dk as 



ki ,fco 



D, 



(lla) 



with k, = {kii,ki2) the coordinates within the plane 
Ekg{ii) — {k I k = kou + fciiv + fci2w} and fco the coor- 
dinate perpendicular to Eko ■ Here, u = v x w and v and 
w are taken to be orthogonal unit vectors. It follows from 
Eq. ( |10c| i that a combination of the symmetries (|7]l and ( |lla[ ) 
imposes the constraint 



gk. 



gki, + fco 



(lib) 



on gk = Alk. For example, in the case of the tetragonal 
point-group with gk given by Eq. (|6a]i, symmetry ( [TT| ) 
is satisfied for the one-parameter family of planes Ek^{i). It 
is important to note that even though Eq. ( |1 Vo\ has the form 
of an inversion-type symmetry, it is different from the mirror 
symmetries imposed by the crystallographic point group, see 
Eq. d5|. We will explain below that the presence of symme- 
try (fm) can give rise to topologically stable arc surface states. 



III. TOPOLOGICAL INVARIANTS IN NODAL 
NONCENTROSYMMETRIC SUPERCONDUCTORS 

To characterize the topological properties of nodal NCSs 
in three dimensions, we introduce an integer topological in- 
variant and two S2 topological numbers. For this purpose it 
is convenient to work in the off-diagonal basis ( fTO] ) and to 
adiabatically deform 'H(k) into a flat-band Hamiltonian. The 
off-diagonal block Dk, which is in general non-Hermitian, 
can through a singular-value decomposition be written as 
Dk ~ ^k^k^k' where J7k and Vk are unitary matrices and 
Sj. is a diagonal matrix with the positive eigenenergies 



Al.k = ^/(e,[)^ + (A^)^ 



A 



2,k 



\/(e 



-^2 + (A^)2 



(12a) 
(12b) 



of H{k) and thus of H{k) on its diagonal, i.e., E^. — 
diag(Ai k, A2,k)- Assuming that Ai k,A2,k 7^ 0, we adia- 
batically deform the spectrum of Ti{k) into flat bands with 
eigenvalues +1 and —1, which amounts to replacing Ek by 



4 



the unit matrix. Hence, the off-diagonal component of the 
flat -band Hamiltonian is given 



where the path C crosses the Fermi surface at k^. It follows 
that the winding number Wc can be simplified to (cf. Ref.il 11) 



q{k) = UlV^ = 



1 



2Ai,kA2,k 



{ 



{Kb^ - KK Iikl }^o = E E ^g" ( ^'^■.?klk=ko, ) (^k°J' (17) 

\ Ik 



(13) 



where = A1.kiA2.k- As a result of TRS, the 2 x 2 unitary 
matrix (/(k) G t/(2) satisfies i(T2<z^(— k) = (7(k)io'2- 



A. Winding number 

In a three-dimensional nodal superconductor, (jf(k) is de- 
fined only for values of k for which Ai k, A2 k 7^ 0, i.e., for 
all k in the BZ except for those on the nodal lines. Hence, the 
three-dimensional winding number which characterizes 
fully gapped systems, is ill-defined for nodal NCSs. However, 
we can use the one-dimensional winding number W n, which 
is defined in terms of a one-dimensional momentum-space 
loop (or line) integral, to characterize the topology of nodal 
NCSs. To that end, we consider ^(k) along a one-dimensional 
loop (or line) in reciprocal space which does not cross gap- 
less regions. That is, we study the mapping U{2) given 
by g(k) € U{2). Since the first homotopy group of U{2) 
IS TTl [U{2)] ^ 'Zf there is an infinite number of homotopy 
classes of mappings from to U{2), which can be labeled 
by the one-dimensional winding number 

= ^ £ dki Tr [q\k)dkM^)] ■ (14) 

The integral is to be evaluated along the path £ parametrized 
by fc;. We observe that Wc is quantized for any closed loop 
£ that does not intersect with nodal lines. If £ encircles a 
line node, Wc determines the topological charge and hence 
the topological stability of the nodal line. Using Eq. ( pj) we 
find 



Wc = ^(f dkidk, [lndet<7(k)] 
2m Jc 

arg (S^ 
[arg (Ck 



1 

2n 
1 

2n 



c 

dki dki 
dki dki 



Al\\^\ 



zA+)+arg (C +*Ak)] 



(15) 



where, in going from the first to the second line, we have used 
the identity |dct(7(k)| = 1. 

Assuming that the energy scales of the gap functions are 
much smaller than those of the normal-state dispersions (i.e, 
for a weak-pairing superconductor), we can rescale with- 
out changing Wc, such that the gaps are nonzero only within a 
small neighborhood of the Fermi surface sheets [cf. Eq. ([12])]. 
In this limit we find that the phase arg(^j^ + iAJ^) is constant 
far away from the Fermi surface of helicity v = ± and that it 
jumps by 



^sgn(9fc,eklk=koJsgn(A;;o; 



where the set of points is given by the intersection of the 
path £ with the Fermi surface for helicity v. Thus Eq. ( [T7| ) 
shows that for a weak-pairing NCS Wc is completely deter- 
mined by the phase structure of the superconducting gaps A^ 
and Aj^ in the vicinity of the positive-helicity and negative- 
helicity Fermi surfaces, respectively. Moreover, it follows 
from Eq. ([TTJi that nodal lines necessarily carry nontrivial 
topological charge Wc 7^ 0, irrespective of the particular form 
of the band structure or the crystal point-group symmetries. 



B. Two-dimensional topological invariant 

The two-dimensional K2 topological invariant introduced 
in this subsection is defined only for NCSs satisfying symme- 
try (fTTb. We therefore consider a Hamiltonian ■H(k) of the 
forrnQ that is invariant under symmetry ([TTJ for the plane 
-Efcg (u). Furthermore, we assume that the nodal lines of 'H(k) 
do not cross Ekg{\\). Hence, ■H(k) restricted to the plane 
Eka (u) can be regarded as describing a two-dimensional fully 
gapped superconductor invariant under both TRS [i.e., sym- 
metry ( [TT| ] and PHS [i.e., the combination of chiral symmetry 
and symmetry ([TT|]. Such a two-dimensional system belongs 
to symmetry class Dili and its topologic al chara cteristics are 
described by the K2 topological invarianP^EIIIEIl 



n 

ki2=Q.,-!r 



Pf [^c^2g(7r, fcj2,fco)] 
Pf [zc72(?(0,fc,2,fco)] 



X e 



(18) 



where g(k) represents a lattice regularization of g(k]^ and 
= —1 (+1) indicates a topologically nontrivial (triv- 
ial) character In Eq. ( [18| l, we have assumed that the coordi- 
nates (fcii, ki2) within the plane Ekg{u.) are chosen such that 
Ki = (0,0), (tTjO), (0,7r), and (7r,7r) are left invariant un- 
der symmetry ( |TT] i. Note that at these points i(j2q(Ki, fco) is 
antisymmetric, so that the Pfaffian is well-defined. 

As before, we consider the weak-pairing limit and set the 
gaps AJ to zero far away from the Fermi surfaces. Using 
Eq. ( fTT] ), we find that in this approximation the exponential 
factor in Eq. ([TSj reduces to 



Yl Yl «sgn(A;;o 5fc,,^l^|k^kO 



(19) 



(16) 



where the set of points S'^ is given by the intersection of 
the Fermi surface for helicity v with the integration paths 
kj : (0,0) — ^ (ti'jO) and ki : (0,7r) — ^ (7r,7r). Assuming that 
the Fermi level does not cross the positive-helicity or negative- 
helicity bands at (Ki,fco), we find that Pt[ia2q(Ki, kg)] is 
either +1 or —1 depending on whether the helicity bands at 
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(K;, ko) are occupied or unoccupied. As a result, the in- 
variant for H(k) simplifies to 



sgn (A+ 



Sgn 



(20) 



where sgn(A e; ) denotes the sign of the gap on the Fermi line 
given by the intersection of Ekg (u) with the positive/negative- 
helicity Fermi surface. Observe that Eq. pO] ) does not depend 
on the lattice regularization and hence is valid also in the con- 
tinuum limit. 

For a NCS that is invariant under symmetry ([TTJ for a one- 
parameter family of planes i?fc(u), with for example e R, it 
is possible to assign a topological charge to the line nodes, 
provided that there are fully gapped regions in momentum 
space separating different line nodes. For concreteness, let 
us consider the situation where a nodal line is located within 
the plane Ek^koi^)- The W12 topological charge of this nodal 
ring can be defined as 



(21) 



where and Ef,- represent two planes that are located 
on either side of the nodal line. For the tetragonal point- 
group [with gk of the form (|6a|], the two-dimensional 
^2 number can be defined for the one-parameter family of 
planes Ek^{z). The nodal rings, which are centered around 
the fc^-axis, consequently carry a nontrivial S2 topological 
chai-ge = iVjP ^^^7V|P , where kp = ^/Zm^I/W 



and a — zt distinguishes between the topological charges of 
the nodal rings in the upper (k^ > 0) and lower (k^ < 0) 
half-spaces. 



C. One-dimensional topological invariant 

Finally, we also introduce a one-dimensional 7^2 invari- 
ant that characterizes the topological properties of 'H(k) re- 
stricted to a time-reversal-invariant loop (or line) £, which is 
mapped onto itself under k — k. Similarly to Eq. ( [T8] l, 
the one-dimensional S2 topological number can be conve- 
niently define d in term s of the lattice-regularized version of 
g(k), Eq. ^ W^m 



N, 



1D_ Pf[»(72g(K2)] _ = 
Pf [lCT2qiK,f 



ii'^dfc,Tr[gt(k)9,,q(k)]^ ^22) 



where we have assumed that C does not cross the nodal lines 
and Ki and K2 denote the two time-reversal-invariant mo- 
menta on the path £. Repeating similar steps as in the previ- 
ous subsection, we find that for a weak-pairing superconduc- 
tor N)P simplifies to 



NjP^ sgn (A+)sgn (A^), 



(23) 



where sgn(A^) represents the sign of the gap at the points 
given by the intersection of C with the positive/negative- 
helicity Fermi surface. 



D. Topological criteria for the existence 
of zero-energy surface states 



As a consequence of a bulk-boundary correspondence,E^E2l 
a nontrivial value of any of the three topological invariants 
( [l4l ), ( [T8] l, and (|22]) signals the appearance of zero-energy 
states at the surface of the NCS. That is, Wc ^ leads to 
surface flat bands, = — 1 gives rise to arc surface states, 
and NjP = — 1 results in Majorana modes at time-reversal- 
invariant momenta of the surface BZ. In the following, we dis- 
cuss in detail the topological criteria for the existence of these 
three types of surface states. For that purpose, we denote the 
coordinates parallel (perpendicular) to a given surface of the 
NCS by r|| (r±) and the corresponding momenta by k|| (fc_L). 

a. Surface flat bands: The appearance of surface flat 
bands can be understood by considering a continuous defor- 
mation of the closed integration path C of Wc, Eq. ( [T4] l, into 
an infinite semicircle, such that the diameter of the semicir- 
cle contains the line (fcj^,k||), with fcj^ e R and k|| fixed. 
This path deformation does not alter the value of the integral, 
as long as no nodal line is crossed while deforming C; if a 
nodal line is crossed Wc changes by ±1. As can be seen from 
Eq. ( [TtI ), the integral along the arc of the semicircle is zero 
and hence the integral along the diameter is given b)I^ 



Wi 



(k|l) = - 



[sgn(Ai:^J-sgn(A;^^ 



(24) 



where \iF.± — (fc_L.±,k||) and k^.i — (A:_L.±,k||) satisfy 
0, respectively, with ]<.f.± (kF.±) cor- 



responding to solutions with positive (negative) signs of the 
Fermi-velocity component perpendicular to the surface. The 
subscript {Imn) in Eq. (|24]i parametrizes the direction per- 
pendicular to the surface. From Eq. (|24]l it follows that zero- 
energy surface states occur whenever Vt^(;mn) (ky ) 7^ 0, which 
corresponds to regions of the surface BZ that are bounded by 
the projections of the nodal rings of the bulk gap. 

b. Arc surface states: Consider an NCS satisfying sym- 
metry ( fTTj i for ^^^(u), with fc e R, and with two nodal rings 
carrying nontrivial Z2 topological charge, Eq. pT| ). Any two- 
dimensional subsystem i?fc(u) lying between these two line 
nodes can be viewed as a time-reversal-invariant topological 
superconductor in class DHL When these two-dimensional 
subsystems are terminated by a boundary, there appear heli- 
cal Majorana edge states. These Majorana modes cross zero 
energy on the surface BZ, somewhere in between the projec- 
tions of the two nodal rings. Thus, there is a one-dimensional 
arc of zero-energy states connecting the projections of the two 
nodal rings. If the projected nodal rings have a finite over- 
lap in the surface BZ, cancellation occurs and no arc surface 
states are expected. By symmetry, the arc surface states al- 
ways appear at those surface momenta that are invariant under 
symmetry ( [TT| i. It is interesting to note that these arc surface 
states are reminiscent of the Fermi arcs that appear at the sur- 
face of Weyl semimetals, which have recently been discussed 
in the context of the A phase of ^He^^' pyrochlore iridates,!^ 
and topological insulator multilayers."^ 
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c. Majomna surface states: We observe that the gap AJ^ 
(with ly — ±.) has the same values at any two momenta related 
by TRS. Thus, for a given time-reversal-invariant surface mo- 
mentum Kii, the one-dimensional W12 number (|23|) read^ 



=sgn(A+ )sgn(AK,.J, 



(25) 



where ii.F,± = (-f^_L±,K||) satisfies 



of the Fermi-velocity component c'ifx I k-K ± 



and the sign 
is assumed 

to be positive. Consequently, iV^^ = — 1 indicates the pres- 
ence of Kramers -degenerate Majorana surface modes at the 
time-reversal-invariant momentum K|| of the surface BZ. 

Before closing this subsection, we note that the topologi- 
cal protection of boundary modes in nodal NCSs is weaker 
than the protection of surface states in strong topological 
insulators or superconductors with a full gap in the bulk. 
Nevertheless, the surface states in nodal NCSs are expected 
to possess a certain robustness against disorder, similar to 
the zero-energy edge states in graphene-^or in da;2_j^2-wave 
superconductors.^^ For example, we expect that weak surface 
roughness will result in a broadening of the surface states de- 
scribed above, but the surface spectral function will still dis- 
play sharp peaks at the bound state energies. Deeper examina- 
tion of the effect of disorder is beyond the scope of this paper, 
and is left to later pubhcations. 



IV. BOUND STATE SPECTRA 

In this section, we use quasiclassical scattering theory to de- 
rive conditions for the existence of surface bound states. We 
will see that the quasiclassical criteria for the appearance of 
zero-energy surface states are in perfect agreement with the 
topological criteria given in the previous section. Hence, this 
provides a verification of the bulk-boundary correspondence 
that relates nontrivial topological characteristics of quasipar- 
ticle wavefunctions in the bulk to the existence of zero-energy 
states at the surface. 



A. General condition for bound states 

The wavevector component ky parallel to an ideal surface 
is a good quantum number due to translational invariance, and 
so we can construct the bound-state wavefunction for each 
point in the surface BZ independently. For a given k||, the 
usual quasiclassical method^* is only applicable if there are at 



most two solutions in each helicity sector. If there are more 
than two solutions in either sector (for example for a concave 
Fermi surface) we are unable to uniquely determine the coef- 
ficients of the different spinors appearing in the bound-state- 
wavefunction ansatz for the given k|| . Our approach is there- 
fore reasonable in the limit of weak to moderate SOC, where 
the concavity of the Fermi surface is only relevant over a small 
fraction of the surface BZ. 

For given k|| we must therefore in general consider four 

wavevectors k± — (fcj^ -|-,k||) and k± ~ (fcj^^±,k||), which 



satisfy 



= 0. We classify a solution as propagating 



if kj_ is £eal and as evanescent otherwise. In the former case, 
k± and k± correspond to solutions with opposite signs of the 
Fermi-velocity component perpendicular to the surface. We 
assume that the NSC is located in the half-space r± > 0. 



1. Propagating solutions on both Fermi surfaces 

If there are propagating solutions on both the positive- 
helicity and the negative-helicity Fermi surfaces we have the 
wavefunction ansatz 



f(k||;r) = ^vI/,(k||;r)e^■ 



I^=± 



*,(k||;r)= J2 a.(k)^.(k)e*^ 



ik±r± „ — K'Cr± 



(26a) 



(26b) 



k=k„,k„ 



where the positive-helicity and negative-helicity spinors are 
given by 



' |ikl+ij 



respectively, with 



-ftSl^k^ 7k+) , (26c) 
1, 7k, SI^O"^' ^26d) 



7k" 



£;-zsgn(4, (k))J|A±P-£;2 



A±p-£;2 



, (26e) 
(26f) 



and Vp j_ (k) is the component of the v — ± helicity Fermi 
velocity normal to the surface. 

A bound state is realized when it is possible to choose the 
coefficients ai/(k) in Eq. (j26| such that the wavefunction van- 
ishes at the surface, i.e., $^||;r)|r^=o = 0. This is equiva- 
lent to the condition 



= (7k_ - \J{il - 7k+J I + I + ^L) - K + ( - iL + *^L) 

+ (7k_ -7t)(7k+, -7£-J - K+^^VilkJ+^U 

X 1 + ^L) ( - It + )-{-it+ ) (Ilk J + ii_ 



(27) 
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Solutions of this equation satisfying \E\ < 
min{|Aj |,|Ai |} are the bound-state energies. Fo- 

cusing on zero-energy solutions that occur within a finite 
region of the surface BZ, we find three possibilities for such 
states: 



(i) sgn(Ak_ A~ J = -1 and sgn(A+^A±^; 

(ii) sgn(Ak_ A~ J = +1 and sgn(A+^A±^; 

(iii) sgn(Ak_ A~ J = -1 and sgn(A+^A±^; 
-1. 



-1, 
-1, 

-1 and 



sgn(Ak_A±^) 



These conditions agree perfectly with the topological criterion 
( p4j l. That is, scenarios (i) and (ii) correspond to topologically 
protected singly degenerate zero-energy states with winding 
number W(imn) = ±1, while (iii) gives doubly degenerate 
states of winding number W^(imri) = i2. We note that a 
careful examination of the terms involving Ik can also yield 
zero-energy dispersing states at isolated points or in a line. 
These states include, but are not limited to, the topologically 
protected Majorana modes and arc surface states introduced 
in Sec. Ill D The general form of this condition is rather cum- 



bersome, but it simplifies significantly in the limit of weak 
SOC and so we defer discussion to this case. 



2. Propagating solutions on only one Fermi surface 

In the case that there are propagating solutions only on the 
negative-helicity Fermi surface, the negative-helicity compo- 
nents of the wavefunction ansatz remain as above, but the 
positive-helicity components are now written 

vl'+(k||;r) = {a+(p)0(p)e-'«p^- 

+ 5+(p)^(p)e''>''-^ } e^P^'-^ , (28a) 

where p is chosen from k+ and k-|_ such that the imaginary 
part of px^ is positive, the spinors 0(p) and 0(p) are defined 
by 



(P) = 



^1 II 



and we have 



1 'pT^"p 'P ' "P p -p 

' iipi+i5' iipi+'r p' f 



r„ = 



A^ 



r„ = 



A^ 



4 \/(A+)^-E^. 



(28b) 
(28c) 

(28d) 
(28e) 
(28f) 



Using the same argument as above, we find the condition 



rp)(7k_ 



(iipi + (iik J + ) - (^p + »^^) ( - + 



(|ipl + ';)(|ii^_l + ^L)-(^P + *'p)(-'L+*'c ) 



(29) 



for the formation of bound states. Unlike Eq. ( |27] l, this condi- 
tion only allows for the existence of singly degenerate zero- 
energy surface states which occur whenever sgn(Aj^ ) ~ 

— sgn(A~ )'P Again, we find that this criterion matches with 

the topological one, Eq. (|24]|. 



3. Limit of weak spin-orbit coupling 

Because of the difficulty of working with spin-orbit split 
Fermi surfaces, many theoretical studies assume that the 
SOC is weak so that the spin splitting of the Fermi surfaces 
can be ignored.'^SEEEl This limit can be directly obtained 
from Eq. ( 27 1 by setting k_|_ = k_ = k and k-|_ = k_ = k to 



yield the compact bound-state condition 

0=(7E+-7k)(7j -7^)(|lk||l£|-lk-l£) 

+ {\ - 7k ) hi - 7k+) (Ilk| llsl + Ik • Ik) . (30) 



This result was previously presented in Ref.TO'and is utilized 
here to obtain the surface bound-state spectra for the O and 
Trf point groups. Although describing a physically idealized 
situation, Eq. ([30| is useful as it significantly simplifies the 
discussion of the zero-energy dispersing states. In particular, 
we find that such states are possible if (i) Ik • Ij^ = |lk| jlj^l ™d 
sgn(A+) = sgn(Ar) =^sgn(A±) =sgn(Ak)or(ii) lk-lk = 



llklll^l and sgn(A+) = sgn(A±) + sgn(A^) 



sgn(Ar). 



Scenario (ii) includes, but is not limited to, the topological cri- 
teria ( |20] i and ( |25| ) which desrcibe arc surface states and the 
Kramers-degenerate Majorana modes, respectively. In par- 
ticular, the antisymmetry of the SOC vector ensures that the 
latter state is realized at the zone center In contrast, states 
satisfying scenario (i) do not fall into either topological cate- 
gory: They cannot be Kramers-degenerate Majorana modes, 
since the requirement that Ik • Ij^ = |lk||lkl is never satisfied 
at the zone center; nor can they be arc surface states protected 
by a %2 number, as the condition on the gap signs implies 
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FIG. 1. (Color online) Comparison of the Fermi surfaces in the x-z 
plane of the approximate dispersion Eq. j32| ) (solid lines) and the 
exact dispersion Eq. 1 31 1 (dashed hnes) for = 0.1. 



that any plane containing k and k must also contain a gap 
node. We note that the condition on the gaps is equivalent to 
sgn(/(k)) = — sgn(/(k)), which is notrealized together with 
the condition on the polarization vector for any of the systems 
considered in this paper, i.e., all examples of zero-energy dis- 
persing states presented below satisfy scenario (ii). 



B. Surface states for the tetragonal point group Cav 

In this section we analyze the bound states of the C^^ point 
group for finite spin-orbit splitting of the Fermi surface. As- 
suming a spherical Fermi surface in the limit of vanishing 
SOC we find the normal-state dispersion at finite A, 




-1.0 -0.5 0.0 03 1.0 



ky I kf 



kyl kf 



no. 2. (Color online) Surface bound-state spectra at the (100) 
face of a Cav point-group NCS as a function of surface momen- 
tum k|| = {ky,kz) with (a) (s -f p)-wave, (b) {d^2_y2 + /)- 
wave, and (c) (d^y + p)-wave pairing symmetry. Here we set 
Afc_F/(2/x) = 0.1 and r = 0.5. The color scale indicates the energy: 
black represents zero energy while yellow represents the maximum 
energy iSmax- The black (gray) line shows the extent of the projected 
negative-helicity (positive-helicity) Fermi surface, (d) Winding num- 
ber W(ioo), Eq. \24\ , at the (100) face corresponding to the same pa- 
rameters as in panel (c). Black (white) indicates M^(ioo) ~ +2 (—2), 
dark blue (gray) corresponds to H^(ioo) = +1 (—1), while light blue 
is VK(ioo) ~ 0. The red dashed (green solid) lines represent the nodal 
lines on the negative-helicity (positive-helicity) Fermi surface. 



4l 
kp 



H±Xjkl 



U2 



(31) 



where kp = ^J^hn[iJW . This implies a concave negative- 
helicity Fermi surface near k^ — ±1^^!, and so our quasi- 
classical method may not be well-defined across the entire 
surface BZ. In order to avoid this problem, we approximate 
the negative-helicity and positive-helicity Fermi surfaces by 
oblate and prolate spheroids, respectively. 



kl - 



(32) 



where 



kp — kp 




(33) 



and they only touch along the line k^ = ky ^ 0. As can be 
seen in Fig.[T] the Fermi surfaces of the approximate and exact 
dispersions agree very well for moderate Xkp/2fi < 0.1. To 
determine the bound-state spectra at given k||, we need only 

find k,y, ky, and the sign of Vpj_(k) at these points so that we 



gives the radius of the helical Fermi surfaces in the kz = 
plane. Our dispersion, Eq. (32i, qualitatively captures the 

salient features of the true dispersion Eq. (31 1: In the Sees. IIID and IV A3 we find that this arc lies along the line 

kz — plane the Fermi surfaces are circles of different radii, ky — 0, i.e., at surface momenta that are left invariant under 



expect that obtaining these parameters from Eq. ( 32 1 should 
give excellent agreement with the exact results across most of 
the surface BZ. 

In Figs. [2|a)-(c) we present the surface bound-state spec- 
tra for the (100) face of a point-group NCS for several 
different pairing symmetries. In the case of (s + p)-wave 
pairing we observe dispersing states in the region bounded 
by the projected positive-helicity Fermi surface and the two 
nodal rings at kz ~ ±0.89fcF [Fig. |2|a)]. Since these nodal 
rings carry nontiivial W12 topological charge, in agreement 
with the discussion in Sec. IIID| there are zero-energy arc 
surface statei^^EoIlSEII connecting the projections of the two 
nodal lings. Moreover, in accordance with the analysis of 
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FIG. 3. (Color online) Surface bound-state spectra at the (101) face 
of a Civ point-group NCS as a function of surface momentum ky 
with (a) (s + p)-wave, (c) {d^2_y2 + /)-wave, and (e) {d^y + p)- 
wave pairing symmetry. Here we set AfcF/(2/i) = 0.1 and r = 0.5. 
The color scale is the same as in Figs.[2ja)~(c). The black (gray) line 
shows the extent of the projected negative-helicity (positive-helicity) 
Fermi surface. Panels (b), (d), and (f) show the winding number 
W(ioi) at the (101) face corresponding to the same parameters as 
in panels (a), (c), and (e), respectively. Dark blue (gray) indicates 
W(ioi) = +1 (—1), while light blue is W(m] ~ 0. The red dashed 
(green solid) lines represent the nodal lines on the negative-helicity 
(positive-helicity) Fermi surface. 



symmetry ([TTJ, and where we have Ik • Ij^ = — jlkHlgl- 
The situation is qualitatively similar for the (dj.2_y2 + /)- 
wave case [Fig. |2|b)], although the extra nodes due to the 
d^2_y2 form-factor /(k) modulate the results of the (s + p)- 
wave case. Note that these extra nodes remove the condi- 
tion for the topological protection of the line of zero-energy 
surface states at ky ~ 0; the zero energy state at the zone 
center nevertheless remains a Kramers-degenerate Majorana 
mode. In contrast, Fig. |2|c) shows that we do not find any 
dispersing zero-energy states for the {dxy + p)-wave pairing, 
but instead there are zero-energy flat bands in several regions 
bounded by the projected line nodes of the positive-helicity 
and negative-helicity gaps. The zero-energy states lying out- 






r=il+2gJ3) r=il+g^/2) r=l 

topological phase topological phase topological phase 

transition, marginal transition, nodal transition, marginal 

nodal points lines touch nodal points 




v=±l ^ 



fully 
gapped 



Wc=±l 




Wc=±l 




v=0 



gapless with 
nodal lines 



gapless with 
nodal lines 



fully j~ 
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FIG. 4. (Color online) Schematic phase diagram for an NCS with 
cubic point group O as a function of the ratio r = As/At of the 
singlet and triplet gaps. Here, the second-order SOC g2 [see Eq. l |6b[ >] 
is assumed to be negative. 



side the projected positive-heUcity Fermi surface (light gray 
line) are the singly degenerate "time-reversal-invariant Majo- 
rana states" found in Ref. TS] and are associated with a non- 
trivial winding number of W(ioo) = ±1 as shown in Fig.[2|d). 
In contrast, the zero-energy surface states lying inside the pro- 
jected positive-helicity Fermi surface are doubly degenerate 
and have winding number W(ioo) = ±2. These states occur in 
the region where the gap has predominantly singlet character, 
and hence are due to the same mechanism as the zero-en ergy 
surface states in a pure dxy-v/ave superconductor.E^Eeissj 

The bound states at the (101) surface shown in Fig.[3]dis- 
play a much more interesting topological character We first 
consider the results for the (s + p)-wave case [Fig. [3|a)], 
which are qualitatively identical to those previously obtained 
in Ref. [To] for vanishing spin-orbit splitting of the Fermi sur- 
faces. Specifically, we find that flat zero-energy bands occur 
within the projected nodes of the negative-helicity Fermi sur- 
face where the sign of the negative-helicity gap reverses be- 
tween the forward- and backward-facing halves of the Fermi 
surface. These zero-energy states are associated with a fi- 
nite winding number W^(ioi) = ±1; the variation of T4^(ioi) 
across the surface BZ shown in panel (b) clearly demonstrates 
the non-zero topological charge associated with the nodal 
rings of the negative-helicity gap. Like for the (100) face, 
the projections of these topologically charged nodal rings 
are connected by arc surface states. The presence of higher 
angular-momentum harmonics [Figs. |3jc)-(f)] results in the 
appearance of additional regions of zero-energy states due to 
the nodes of both the positive-helicity and negative-helicity 
gaps. All of these states correspond to a winding number of 
W^(ioi) = ±1, as can be seen by comparing the bound state 
spectra. Fig. [3]^c) and (e), with the winding number calcula- 
tions. Fig. [3jd) and (f), respectively. 
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FIG. 5. (Color online) (a)-(e) Surface bound-state spectra for the (100) face of a NCS with point group O, g2 ~ —1.5, and A = 0, as a 
function of surface momentum k|| for (a) r = 0, (b) r = 0.15, (c) r = Vc = 0.25, (d) r = 0.6, and (e) r — 0.9. The color scale is the same as 
in Figs.[2ja)-(c). (f)~(j) Same as panels (a)-(e) but for the (110) face, (k)-(o) Same as panels (a)~(e) but for the (111) face, (p)-(t) Winding 
number for the (111) face corresponding to the same parameters as in panels (k)-(o). Dark blue (gray) indicates M^(iii) = +1 (— 1), 

while light blue is W(iii) = 0. The red dashed lines represent the projections of nodal lines. 



C. Surface states for the cubic point group O 



In the cubic point group O there are pronounced changes in 
the nodal structure of the negative-helicity gap as the singlet- 
to-triplet ratio r = As /At is varied, which are reflected in the 
surface bound-state spectrum. To simplify the discussion, we 
assume a spherical Fermi surface, negligible spin-orbit split- 
ting, and finite 52 < O.'^^ A schematic topological phase dia- 
gram of this NCS is presented in Fig. |4] For r > 1 the sys- 
tem is fully gapped and topologically trivial. Reducing r, we 
find that at r = 1 point nodes appear in the negative-helicity 
gap at k = fcf (1, 0, 0) and equivalent points. This is a Lif- 
shitz transition^ at which the topology of the Bogoliubov- 
de Gennes quasiparticle spectrum changes, but the symme- 
try of the ground state remains unaltered.*' Upon further 
lowering r, these point nodes develop into nodal rings with 



non-trivial topological charge. A critical value is reached at 
r = Tc = 1 + (72/2, where the nodal rings touch each other 
and reconnect in a different manner, i.e., there is a change in 
the topology of the nodal structure itself, with the rings now 
centered about k = {kp/^/3){l, 1, 1) and equivalent points. 
Finally, at r = 1 + 2(72/3 there is another Lifshitz phase tran- 
sition to a fully gapped phase with topologically non-trivial 
order characterized by the three-dimensional winding num- 
ber I' ~ il.** The nodal structure at this transition point is 
marginal in the sense that it is topologically trivial and not pro- 
tected against decay either into the fully gapped topologically 
non-trivial state or into the topologically stable gapless phase, 
i.e. point nodes do not posses any topological protection in a 
three-dimensional NCS. Experimental signatures of the topo- 
logical phase transition will be discussed in Sec. |V B| In the 
following discussion of the bound states we take 32 — —1.5, 
which implies line nodes for < r < 1, point nodes for 
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FIG. 6. (Color online) (a)-(e) Surface bound-state spectra for the (100) face of a NCS with point group Td and A = as a function of the 
surface momentum ky for (a) r = 0, (b) r = 0.2, (c) r = = 4\/2/9, (d) r — 0.75, and (e) r — 0.9. The color scale is the same as in 
Figs. [2ja)-(c). (f)~(j) Same as panels (a)-(e) but for the (110) face, (k)-(o) Same as panels (a)-(e) but for the (111) face. (p)~(t) Winding 
number for the (111) face corresponding to the same parameters as in panels (k)-(o). Dark blue (gray) indicates lV(iii) = +1 (— 1), 

while light blue is W(iii) = 0. The red dashed lines represent the projections of nodal lines. 



r — 0,1, and a fully gapped state for r > 1. Surface bound 
states only occur for r < 1. 

We first consider the bound states at the (100) face 
[Figs. |5]^a)-(e)], which were partially examined in Ref. 27 
for r > Tc- Starting at r = we find bound states for all 
|kj| I < kp except at the projected point nodes. Of particular 
note are the dispersing zero energy states at (fcij[ , fc2.|[ ) = 
(0,0), (A:i?/\/3, 0) and equivalent points, where the condi- 
tion Ik 



-Ij^ is satisfied. 



Ill 



however, we de- 



From Sec. 

duce that only the zero-energy state atlfie zone center is a 
topologically protected Majorana fermion. The other four 
zero-energy states, in contrast, only form for g2 < —1 and 
hence are not topologically stable. The apparent zeros close to 
{kp/VS, kp/VS) and symmetry-related points in (a) are the 
projections of the point nodes. Upon increasing r, the nodal 
rings on the negative-helicity Fermi surface lead to regions 



within the projected Fermi surface where bound states do not 
occur, i.e., the white space in Fig. |5jb). Note that we have 
W^(ioo) = everywhere since the integration path £ passes 
through none or two nodal rings, the contributions of which 
cancel. The reconnection of the nodal rings at r — Tc splits 
the bound states into disconnected regions at the zone center 
and near the edges of the projected Fermi surface. Note that 
the Majorana fermion mode at the zone center survives up to 
r = 1. As shown in Figs.|5];f)-(j), the situation for the (110) 
face is essentially similar, although here the Majorana fermion 
state at the zone center only survives up to the topological Lif- 
shitz transition at r j-^- 

The bound states at the (111) face are qualitatively dif- 
ferent due to the presence of topologically protected zero- 
energy bands. We show the bound state spectra and the as- 
sociated maps of the winding number W^(iii) in Figs. [5|k)- 
(o) and Figs. |5]^p)-(t), respectively. Unlike for the other two 
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FIG. 7. (Color online) Tunneling conductance spectra for the (100) and (101) interfaces of a CUv NCS with r — 0.5, Z — 3, and (a) (s + p)- 
wave gap symmetry [/(k) — 1], (b) {d^2_y2 + /)-wave gap symmetry [/(k) = (fc^ — ky)/k%], and (c) {dxy + p)-wave gap symmetry 
[/(k) = 2kxky/kp]. In all panels the thick lines show the results for finite spin-orbit splitting Xkp/2n = 0.1, while the thin lines show the 
results for degenerate helical Fermi surfaces. 



surfaces, we find that at r = the energy of the dispersing 
states reaches the gap so that we have regions within the pro- 
jected Fermi surface where bound states do not occur. For 
finite but small r, the dispersing bound states break up into 
six lobes, with zero-energy flat surface bands within the pro- 
jected nodal lines close to the center and edges. Like for the 
point group, the zero-energy states are associated with a 
winding number — ±1. Increasing r we see that the 

zero-energy surface states grow in extent while the dispersing 
states shrink, until at r Tc the nodes in the gap touch and 
the zero-energy states suiTound the dispersing states. The dis- 
persing states vanish at r closer to 1 when the projected nodes 
no longer overlap, and we hence find only zero-energy states 
inFig.jSjo). 



D. Surface states for the tetrahedral point group Ta 

The tetrahedral point group Td also shows topological 
phase transitions in the nodal structure upon varying r, which 
we again discuss for a spherical Fermi surface and vanish- 
ing spin-orbit splitting. At r > 1 the system is fully gapped 
and topologically trivial. Upon reducing the singlet-to-triplet 
ratio to ?■ = 1, a Lifshitz transition occurs with the appear- 
ance of twelve point nodes in the negative-helicity gap at 
k = {kp/^/2){l, 1, 0) and equivalent points; further reducing 
r we find topologically charged nodal rings centered about 
these points. These rings grow as r is lowered, eventually 
touching at r — Tc — 4v^/9. For r < rc the nodal rings re- 
connect in a different way so that they are now centered about 
the fourteen zeros of Ik at k = kpil, 0, 0), (fcF/\/3)(l, 1, 1) 
and equivalent points. In the limit of a purely triplet gap we 
find point nodes at these locations. 

The bound states at the (100) surface are shown 
in Figs. |6|a)-(e). Unlike for the cubic point group O, there 
are no zero-energy dispersing states for this surface; the ap- 
parent zero-energy states at r ~ in Fig. [6|a) are in fact 
nodes. For r > these points become lines, clearly visible 
as the boundaries of the white space. As r is increased past 
the dispersing bound states separate into disconnected regions 



of the BZ, which vanish at r = 1. We now turn to the bound 
states for the (110) surface, which are presented in Figs.|6|f)- 
(j). At r < we find lines of dispersing zero-energy states 
which connect the projected point nodes (for r = 0) and nodal 
rings (for r > 0) of the negative-helicity gap. These arc states 
are not topologically protected by the mechanism discussed in 



Sec. Ill D however, since for this point-group symmetry the 
condition Eq. ( 1 lb i only holds for planes which intersect line 



nodes of the gap. We nevertheless note that the zero-energy 
state at the surface BZ center is a topologically protected Ma- 
jorana fermion which is present for all r < 1. 

Like for the cubic point group O, topologically protected 
zero-energy bound states are found at the (111) surface for 
< r < 1. The bound state spectra are plotted in Figs.|6|k)- 
(o), with the variation of the winding number shown 
in Figs.|6|p)-(t). The main difference to the case of the point 
group O is the much more limited extent of the dispersing 
bound states, which are almost entirely absent for r > rc- 
At small r we also observe arc surface states connecting the 
topologically charged nodal rings near the edge of the pro- 
jected Fermi surface, but again these states are not topologi- 
cally protected by the mechanism discussed in Sec.|IIID] 



V. TUNNELING CONDUCTANCE AND TOPOLOGICAL 
PHASE TRANSITIONS 



In this section we discuss the calculation of the tunneling- 
conductance spectra for a normal-metal-NCS junction. Such 
tunneling experiments are an important test for the pairing 
symmetry of unconventional superconductors. In paiticu- 
lar, they can confirm the existence of the zero-energy sur- 
face flat bands, which are evidenced by a sharp peak in the 
low-temperature conductance at zero 

bias.UaOi We also dis- 
cuss possible signatures of the topological phase transitions. 
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A. Tunneling conductance 

At zero temperature the tunneling conductance as{eV) for 
tunneling from a normal metal without SOC into an NCS 
under bias voltage V is g iven by t he generalized Blonder- 
Tinkham-Klapwijk formulcPSHEl 

a5(ey) = ^^{l + |a.,.,(kji) 



l&...'(k||)P} 



(34) 

where aa-.a-'0^\\) and 6o-.cr'(k||) are the spin-resolved Andreev 
and normal reflection coefficients for electron injection into 
the NCS, respectively. These coefficients appear in the wave- 
function ansatz describing the electron-injection process, 

*,(k||,r) =7/^<(k||,r)e(-ri) + ^>(k||,r)e(ri). (35) 

In the normal metal the wavefunction is written as 



ik-r 



+6^,<,'(kj|)V'e,a'e'''"" 



(36) 



with the electron and hole spinors ■0e CT = |(1+(t, 1— cr, 0, O)"'" 
and -0/1. o- = i(0, 0, 1 + fj, 1 — a)'^, respectively. In the case 
that there are propagating solutions in both helicity sectors, 
the wavefunction in the NCS is written as 

V'>(k||,r) = [c„,„V»(k„)e'''"-'' + d^,„Vn(k„)e'''" 

n=± 

(37) 

where the spinors are defined as in Sec. IV A 1 if there 
are propagating solutions only in the negative-helicity sec- 



tor, the positive-helicity components of Eq. ( |37| l are replaced 
by Eq. 



P8a| i. To simplify the discussion we assume that the 
normal metal has a spherical Fermi surface with the same 
chemical potential and effective mass as the NCS. We restrict 
ourselves to the limit where the bias energy is negligible com- 
pared to the Fermi energy so that the energy dependence of 
the wavevectors can be ignored. 

The reflection coefficients are determined by application of 
the interface boundary conditions. The first condition requires 
continuity of the wavefunction across the interface, i.e.. 



rj_=0+ 



(38) 



The second condition enforces conservation of probability 
across the barrier,'^^'^ which is modeled as a (5-function of 
height U. We have 

*sV^>(k||,r)|,^^o+ = {dr^+2ZkF)yj<m,r)l^^o- , 

(39) 

where Z ~ mU/fi^kp is a dimensionless parameter charac- 
terizing the transparency of the interface, and vs is propor- 
tional to the r± component of the velocity operator in the nor- 
mal state of the NCS. For a 6*4^, point-group NCS orientated 
such that the b axis is parallel to the interface and the c axis 
makes an angle a to the interface, we write 



vs = (To (8) (7o dr^ +id-o< 



mX 

' 0-2 -r;^ COS a . 




FIG. 8. (Color online) Variation of the zero-bias conductance-peak 
height at the (111) face as a function of singlet-to-triplet ratio r for 
the point group O. The black, downward-pointing arrows indicate 
topological phase transitions between fully gapped and nodal phases. 
The red, upward-pointing arrow marks an inflection point in os (eV), 
where the character of the nodes on changes qualitatively. In this 
figure we set 32 = — 1, = 3, and T = K. 



In Fig.|7]we present conductance spectra for tunneling into 
a C4i, point-group NCS through (100) and (101) interfaces 
for both vanishing and finite SOC A. We find the tunneling- 
conductance results to be quite robust against the finite spin- 
orbit splitting of the Fermi surfaces. Several features of these 
spectra are noteworthy. For the (100) surface we observe a 
broad hump-like feature in the tunneling conductance for the 
(s + p)-wave and (d^2_y2 + /)-wave pairing states, which is 
a signature of the arc surface states. In the (dxy + p)-wave 
case, in contrast, we find a zero-bias conductance peak well- 
separated from the bulk density of states. In the limit of degen- 
erate helical Fermi surfaces this is due to the doubly degener- 
ate zero-energy states for k|| close to zthpe^ [see Fig. |2|d)], 
while for non-zero SOC the singly degenerate zero-energy 
states near k|| = ztkpBy also contribute. For the (101) surface 
all pairing symmetries show a zero-bias conductance peak, 
which is a key experimental signature of the topologically pro- 
tected zero-energy states.'" '-^ For the cases of (s + p)-wave 
and {d^2_y2 + /)-wave pairing, however, we note that this is 
superimposed on a much diminished hump-like feature. This 
signals the continued existence of arc surface states in these 
systems, in agreement with Fig. [5] 



B. Evidence for topological phase transitions 



(40) 



In Sec. IV C we have demonstrated the existence of topo- 
logical phase transitions as a function of the singlet-triplet ra- 
tio r for an NCS with point group O and SOC vector Ij. given 
by Eq. ( |6b] i with 52 7^ 0. This discussion is physically rel- 
evant for Li2Pdj;Pt3_j;B, for which the SOC strength can be 
tuned by substituting Pd for Pt.^^ The magnitude of the SOC 
interaction in these compounds seems to be directly related to 
r,-'* which suggests that it might be possible to observe topo- 
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logical phase transitions between a fully gapped and a gapless 
phase, or between two gapless phases, in Li2Pd2;Pt3_2.B as a 
function of Pt concentration. 

The most direct way to detect the topological phase tran- 
sitions in this system requires measurements sensitive to the 
low-energy bulk density of states p{Li>). At the r = 1 and 
r = (1 + 2(72/3) Lifshitz transitions, this quantity changes 
from p{uj) = 0, characteristic of a full gap, to the linear de- 
pendence p{lj) (X uj associated with line nodes. Signatures of 
the topological phase transition can also be seen in the con- 
ductance spectra, in particular the zero-bias conductance peak 
at the (111) surface. Indeed, as shown in Fig. [8] the zero-bias 
conductance shows abrupt changes at the boundaries of the 
nodal region. Furthermore, at the critical where the topo- 
logical structure of the nodes changes, there is a kink anomaly 
which is marked by the red, upward-pointing arrow in Fig. [8] 
Tunneling experiments could therefore in principle be used to 
evidence a topological phase transition in Li2Pda;Pt3_2;B. 

VI. CONCLUSIONS AND OUTLOOK 

We have performed a detailed analysis of the topological 
properties of nodal NCSs in three dimensions. Using topolog- 
ical arguments, we have derived general criteria for the exis- 
tence of Andreev bound states at the surface of nodal NCSs. 
Three different types of topologically protected surface states 
have been identified, namely Kramers -degenerate Majorana 
modes, arc surface states, and surface flat bands, whose sta- 
bility is guaranteed by the conservation of the one- and two- 
dimensional ^2 invariants [Eq. ( |22] i and (fTS])], and the wind- 
ing number ( [T4] i, respectively. 

We have independently derived the general criteria for zero- 
energy surface states using the quasiclassical scattering the- 
ory. Furthermore, we have applied this technique to study a 
number of physically relevant manifestations of NCSs. For 



a C4t, point-group symmetry, we have calculated the surface 
bound-state spectra and the tunneling conductance for (s 
wave, ((i^2_y2 + /)-wave, and (d^y -V p)-wave pairing and 
finite spin-orbit splitting. We have shown that the surface 
bound states are in perfect agreement with the variation of the 
topological winding number across the surface BZ. We have 
shown how the surface flat bands manifest themselves as a 
zero-bias conductance peak, while the arc surface states lead 
to a broad, hump-like feature centered around zero bias in the 
conductance spectra. Both features exhibit a pronounced de- 
pendence on surface orientation, which provides characteris- 
tic fingerprints of the topological properties of the system. We 
have also examined the bound-state spectra in NCSs with O 
and Td point-group symmetry in the limit of weak spin-orbit 
splitting and discussed the occurrence of topological Lifshitz 
transitions in these spectra. 

On symmetry grounds the zero-energy surface flat bands 
are expected to appear in any unconventional nodal super- 
conductor preserving TRS. In particular, this should be the 
case for the class-Cl topological superconductors-l^^lES Fur- 
thermore, flat bands are also expected to occur at the surfac e of 
nodal topological superconductors in symmetry class 
(i.e., topological superconductors with TRS that are invariant 
under rotations about one fixed axis in spin space). Besides 
the tunneling conductance, the surface flat bands and arc sur- 
face states also profoundly affect other surface and interface 
properties of NCSs, such as Josephson tunneling,^' the non- 
linear Meissner effect, and surface thermal transport. The in- 
vestigation of these interesting boundary properties are left for 
future work. 
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